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Abstract 



We analyze the one loop corrections to leptonic W and Z decays in an R- 
parity violating extension to the Minimal Supersymmetric Standard Model 
\^ , (MSSM). We find that lepton universality violation in the Z line-shape vari- 

l/^ \ ables alone would strengthen the bounds on the magnitudes of the A' cou- 

f^ ' plings, but a global fit on all data leaves the bounds virtually unchanged at 

Q ' 1-^33^1 — 0-42 and |A23^| < 0.50 at the 2a level. Bounds from W decays are 

less stringent: lAgg^^l < 2.4 at 2a, as a consequence of the weaker Fermilab 
experimental bounds on lepton universality violation in W decays. We also 
point out the potential of constraining R-parity violating couplings from the 



Q^' measurement of the T invisible width. 
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I. INTRODUCTION 

The assumption of R-parity conservation in supersymmetric model-building has long 
been an economical means of (1) avoiding certain phenomenological problems in SUSY 
models [e.g. proton decay), (2) ensuring that the lightest supersymmetric particle is avail- 
able as a cure for the dark matter problem, and (3) reducing the SUSY model parameter 
space. (For recent reviews, see Ref. [Q.) However, the recent discovery of neutrino mass 
at Super-Kamiokande provides improved motivation for R-parity violating extensions 
to the Minimal Supersymmetric Standard Model (MSSM). Detailed analyses of the phe- 
nomenological constraints on such models are thus warranted to quantify the amount of 
R-parity violation permitted by current experimental data. 

In this paper we consider the effects of R-parity violating extensions to the MSSM on 
lepton universality in W and Z decays. The R-conserving sector of the MSSM generates 
lepton universality violations proportional either to the lepton Yukawa couplings (due to 
Higgs interactions) or to the mass splittings of the sleptons (due to gauge interactions). 
Effects due to the Higgs sector will be considered in a future work ||^ but will in general be 
negligible unless tan/3 is quite large 0]. Effects due to gauge interactions are negligible if the 
slepton mass splittings are small. This is the case, for example, in supergravity (SUGRA) 
models with universal soft-breaking scalar masses at the SUGRA scale, in which the mass 
degeneracy is broken only by renormalization group running effects involving small Yukawa 
couplings. In R-parity violating models, however, R-parity violating interactions provide 
additional sources of lepton universality violation which may be significantly larger than 
these smaller effects. 

The R-parity violating superpotential has the following form:[| 

1 ... ^ . . . I „ -■ -■ -^ 

^fi ^ T^^ijkLiLjEk + Xiji^LiQjDk + -\jiJJiDjDk , (1.1) 

where Lj, Ei, Qi, Ui, and Z)j are the MSSM superfields defined in the usual fashion 0, and 
the subscript i = 1,2,3 is the generation index. Since a priori the interactions described 
by this superpotential have an arbitrary flavor structure, we generically expect that in the 
context of W and Z decays these will give rise to lepton universality violations. In this 
paper, we estimate the size of this violation and derive constraints on the R-parity violating 
couplings from LEP and Fermilab measurements of the leptonic observables in W and Z 
decays. 

It is clear that the purely baryonic operator JJiDjDk is irrelevant to our discussion. The 
other two operators may affect W and Z decays at one loop through vertex corrections to the 
WClt^ and ZdidL vertices, with superparticles running in the loop. However, the couplings 
Xijk are already tightly constrained to be at most (9(10~^), as the operator LiLjEk violates 
lepton universality in lepton decays at tree level ^. The constraints on A^ ;, are much less 
stringent. Previous limits cite upper bounds on X'^.^ as large as the gauge couplings {i.e. 



^Because we are neglecting the soft breaking terms, we can rotate the bilinear terms away |||]. For 
possible effects of the soft breaking terms on W and Z decays see Ref. P]. 



as large as 0.5) with the SUSY scale at 100 GeV |T[]. One may thus expect significant 
radiative corrections induced by these couplings. Henceforth we will focus on the effects of 
the operator LiQjDk only. 

It is important to note that very strict constraints on the products of liz/ferent R-violating 
couplings already exist from fiavor- changing processes, e.g. yU ^ 67 constrains lA'^j^A'ajjl < 
4.6 X 10~^ 0. However, these constraints can be easily satisfied by requiring only one of these 
couplings to be very small leaving the other coupling ill-constrained. The W and Z decay 
processes we consider here are fiavor-conserving and involve the same R-violating coupling 
squared. Thus we can constrain the individual couplings rather than their products. 

We emphasize that, in the absence of a complete calculation in the full theory, focussing 
attention on the violation of lepton universality provides a clear advantage over studying the 
effects of R-breaking couplings on the individual lepton-gauge boson couplings separately. 
This is because the R-conserving sector induces significant universal corrections to the lepton 
couplings which depend strongly on the choice of SUSY parameters. These corrections (along 
with the corrections to the hadronic partial widths) cancel when considering violations 
of lepton universality. Thus, the study of lepton universality violation lets us isolate the 
effects of R-breaking interactions without ad hoc assumptions about corrections from the 
R-conserving sector. 

In the following calculations we neglect left-right squark mixing. Left-right squark 
mixing could be large only for the stop. However, since diagrams involving the stop contain 
down quarks with negligible mass, it will be seen that the contributions from these diagrams 
are numerically small (subleading in an expansion in rn^ or m|). Further, due to the chiral 
structure of the R-breaking interactions, two left-right mass insertions would be required 
in the diagram. Thus, such contributions would be further suppressed as long as the mixing 
parameter is perturbatively small. 

Radiative corrections to individual Z ^* II partial widths due to R-breaking interactions 
have previously been considered in Ref. [^ but lacked a consistent treatment of the R- 
conserving corrections. In this paper, we study the violation of lepton universality in W 
and Z decay to isolate the effects of R-breaking couplings and constrain their sizes. In 
determining the limits on R-breaking from Z decay, we perform a global fit to all the 
relevant LEP and SLD observables in which the corrections from both R-breaking and R- 
conserving interactions are parametrized and fit to the data. This provides a consistent 
accounting of R-conserving effects and allows us to improve the existing bounds on the 
R-breaking A' couplings. A companion study of constraints on A' and A" couplings from 



LEP/SLD hadronic observables has been performed in Ref. |TT 



II. LEPTONIC W DECAYS 

The relevant R-parity violating interactions expressed in terms of the component fields 
take the form 



A£^ - \jk 



i^iLdkRdjL + djid^RViL + dlj^u'^idjL 



-{eiLdkRUjL + UjLdkRCiL + dlj^el^UjL) 



h.c. (2.1) 



The one loop diagrams contributing to the decay W -^ eni^i'L are shown in Figs. |A2| and 



"KQ, . At one loop, the neutrino flavor may differ from that of the tree level vertex {i ^ i') 
as a result of the R-parity violating interactions. Since neutrino flavor is indistinguishable 
in the detector, we should in principle sum over all three generations of antineutrino in the 
final state. However, since this is a one loop effect which does not interfere with the tree 
level flavor conserving decay, we will neglect it in our analysis and set i = i' . 
The amplitude of each diagram in Figs. |A2| and |A2| are: 



-NnW. 



C \^ijk\ 



-^4^W^{p + q)eiL{ph^iyiL{q) 
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'1 ^^UjL 



5i 0;0,m^ 



■JL 






' O-kR 



(2.2) 



5i(0;0, 

The expression in the square brackets is the tree level amplitude. The definitions of the 
integrals Bi, C23, and (^24 are presented in the Appendix. In the above expressions Nc = 3 
is the number of colors, and i, j, k are family indices; the final result must be summed over 
j and k to obtain the full correction for final state flavor i. We have set all the down type 
quark masses to zero. The up type quark masses m„, will also be set to zero except for the 
top quark (j = 3, m^g = rrit). 

Combining the expressions in Eq. (|2.2D, with appropriate factors of | for the wavefunction 
renormalizations, we obtain the one-loop shift of the WenViL coupling due to the \[jj^ 
interaction: 



^9: 



ijk 



Sglfk + ^aijl 



S9^ 



I |2 



^9 



(u) 
ijk 



-NclKjk 



-Nc\X[,k\' 



2(724 (0,ms^,„m,-^.J + ^B, (o,m^^,,) + ^B, (o 



' d.jL 



+U (0 



' dkR 



+ -Bi[mu^,m^^^ 



(2.3) 



We have suppressed the external momentum dependence of the B and C functions to sim- 
plify our expressions. The contributions of diagrams involving the down-type quark have 
been combined in Sgljl and those involving the up-type quark in Sgljl- The 1/e poles of 
dimensional regularization cancel separately in each of these combinations so they are finite. 
In the following, we evaluate the size of these shifts for a common squark mass of rrig = 
100 GeV. To facilitate estimations for different squark masses, we provide approximate 
formulae. 



First, we evaluate 5gl^l, the contribution from diagrams involving the top quark ^3. 
For m„3 = rrit = 175 GeV, ruw = 81 GeV, and rriq = 100 GeV, we find: 

^^ = -1.02%|A' P 
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'iSfcl 



(2.4) 



An approximate expression can be obtained by expanding the full expression of 5gl^l 
in powers of m^: 



^9i3k 



N, 



c ,,/ ,2 



(A^\2\'^'i^k 



(47r) 



A 



X 



4(1 -x)' 



+ 



m 



w 



mi 3(1 — x) 



-| X — 1 — (2 — x) Inx !• 
Al — X + hix\ 



(2.5) 



where x = m'^/m'^. For rriq = 100 GeV, this expression is equal to (—1.11 + 
0.09)% lA^g^jp. Compared to the exact result above, we see that the leading order 
approximation is already fairly accurate. 

The contributions of the diagrams with massless quarks are numerically smaller and 
vanish as mw -^ 0. The correction with massless up-type quarks (j = 1,2) is: 



^^ = 0.22% I A', 12 
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^ijk I 



(J = l,2) 



The approximate form to leading order in m'^/m'^ is 



^9^ _ Nc ^y |2 ml, f^ 


31n"^'^ 
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(47r)2' '3''' 9ml \ 


m? 



(2.6) 



(2.7) 



For rrig = 100 GeV, this gives 0.31% |A^ ;,p which suffices for our purpose. 
The diagrams with down-type quarks contribute: 

Jl2li = 0,07%\X' 12 



9 



^ijk I 



(2.8) 



the approximate expression being 

S9^L . . Nc 



Hjk 

9 



(47r 



A 



' |2 



m, 



w 



1 2 ^ ijk 



18 m?^ 



(2.9) 



For m,g = 100 GeV, this gives 0.07% |A^ ;,p to the accuracy shown 



Note that each contribution, Eqns. ( p.5|) , (p7|), and ( pl9|) separately decouples in the limit 
m? ^ oo as they should. Collecting everything together, the shift in the coupling of the 
z-th generation lepton to the W is given by: 
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E 



^aU ^ ^aijl 



-0.95%^ |A:3,P + 0.29%j:\K2k\' + 0.29%j:\Kik\' (2.10) 



where we have summed over all possible generation indices j and k. 



The current bound on lepton universality violation in leptonic W decays from D0 is |13 



— = 1.004 ± 0.019(stat.) ± 0.026{syst.). 

9e 



This places a constraint on 



6 



^93 Sgi 

9 9 ' 



(2.11) 



Note that R-conserving corrections do not contribute since they cancel in the ratio Qt/ 9e- 
The constraint on the R-breaking couplings is thus 



-0.4±1.9(stat.)±2.6(s|/st.) 

= {E|AmP-0.3ElA'mP-0.3ElA^i.P} 

-{EIAmP-o.3EiA;2,P-o.3EiA;,,r}. 



(2.12) 



Of the couplings X'^j^ appearing in this expression, the i = 1 couplings are already well 
constrained from neutrino-less double beta decay, neutrino masses, atomic parity violation, 
and low energy charged current universality. Using the 2cr limits for the individual couplings 
charted in Ref. 113, we find for rripj = 100 GeV: 



eia; 



iifci 



E l-^12fcl 
k 



< 0.00088, 



< 0.0055, 



E l-^13fcl 
k 



< 0.079. 
The i = 3, j = 1 couplings are also well constrained from Rj 



r(r — > 7rz/T-)/r(7r 



Again, using the 2a limits cited in Ref. [|I^ we find 

^|A:,i,P< 0.036 



(2.13) 
(2.14) 



Therefore, we can neglect the A'^^^i, and Ag^^, terms in Eq. ( |2.12|) and obtain 

ElAmr-0.3ElAmr = -0.4±3.2 



(2.15) 
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where the systematic and statistical errors have been added in quadrature. If we neglect the 
32k term with a smaller numerical coefficient, this places a la {2a) upper bound on the 33k 
term: 



^|A;,3,|2<2.8(6.0), 

k 



(2.16) 



which in turn translates into the limit 



|A;,3,|<1.7(2.4). 



(2.17) 



Non-zero values of X'^2k ^i^l weaken this bound. 



As we will see later, Z decay data places a constraint on J2k l-^ssfcP ^^ ^^e ±0.1 level. 
Therefore, for the W decay data to be competitive with the Z decay data, the error must 
be improved by more than an order of magnitude. While the Tevatron Run II may provide 
enough data to improve the statistical error considerably, improving the systematic error 



may prove a challenge |15 



III. LEPTONIC Z DECAYS 

We next consider the effect of R-parity violating interactions on flavor-conserving lep- 
tonic Z decays, Z -^ eneii. Note that the A' interaction in Eq. ( |1.1| ) involves only the 
left-handed lepton field. Therefore, at the one-loop level the right-handed coupling is un- 
affected. Neglecting all down-type quark masses, the amplitudes of the diagrams shown in 
Figs. |A2| and |A2| are 



where 



-Nc\K '^ 



Hjfc 
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Z^{p + q) eiL{p)l^,eiL{q) 



(3a) 
(36) 
(3c) 



cos 6w 
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+2hd^C2i (O, 0, ml; mu,,m 



jLl'^UjL 



m: 



-K 



{d - 2)^24 (O, 



> ^z; ^d. 



'^M, ; ^^ti 



mzC23 0,0,m^;m.- ^,m„,,m„ 



dfcfl' 



{3d) : +hd^ [{d - 2)^24 (O, 



,m^;mji^.^,0,0 



-mlC23 (O, 0, m|; m^^.^, 0, 



m^\m- 



,, iioz, 



(3e) : K^mlCo{o,0 

(4a) + (46) : 2K^B, {0;0,mu^^ 

(4c) + (4ci) : 2/i,^5i(0;m„^.,m4j 



,mu.,m. 



(3.1) 



hf^ = h-Qf sin^ 6 



w, 



hfH = -Qf sin^ 



w- 



(3.2) 



The tree level amplitude is /ig^ times the expression in the square brackets. These corrections 
can be expressed as a shift in the coupling he^: 



5h 



ijk 



shU + sh^ 



ijh 



^^ijk = -^c\\jk\ 



2/i«iC24 0,mii m^ 



^Kjk = -NcWjk\ 



+hd^ {{d - 2)C'24 (m^,,, 0, O) - mlC23 (m^,-„ 0, o)} 

2/1^^(724 (m„^,m4^,m4^) 
-huL {{d- 2)C'24 (mj^^, m„^,, m„J - ml^ss ("^4^, m„., m„J } 



-/iei5ifm„^,,mj^^ 



(3.3) 



Again, the dependence on the external momenta has been suppressed. The corrections which 
depend on the down-type quark have been grouped together in Sh\jl and those that depend 

on the up-type quark in Shi^l- These combinations are separately finite. For simplicity, we 
again evaluate these shifts for a common squark mass of rriq = 100 GeV. 

• We begin with the top quark dependent contribution. For rrit = 175 GeV, rriz = 
92 GeV, and rriq = 100 GeV, we find 



6h) 



(u) 



i3k 



0.63%|A^3fc|l 



(3.4) 



This is well approximated by the leading mz = piece of the expansion in the Z mass: 

Nr 



SK 



(u) 
ijk 



'C ix/ ,2 



2(47r)2 



K,kvn^) 



(3.5) 



where 



^/ N X f 1 , 

Fix) = 1 H In a; 



1 — a; V 1 — X 



(3.6) 



and X = m^/m?-. For niq = 100 GeV, this gives 0.65% jA^g^^p. The subleading terms 



An) 



from the individual diagrams contributing to 5h\^l. are: 



m 

(3c) 
(3e) 



I \/ \-A ■■''Z I- 



(47r)2 

Nc Ki^ , w \2^Z 
fA \9 \ i3k\ 2 

(47r)^ mf 



2m| '' \x 
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fi^) 



(47r)2 '^^3/cl ^2 9K^) 



(3.7) 



where 



fix) 

9{x) 
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1 
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(1 
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X)4 


1 




1 





12(1 



x] 



2x^ — 9x^ + ISo;^ — llx — 6a; In X 
x^ — 6a;^ + 3x^ + 2x + 6x^ In x 



(3.8) 
(3.9) 



The total subleading contribution for nig = 100 GeV is —0.03% |A^j;,p. 

In Ref. [|I0[, the leading and subleading contributions of diagrams (3b) and (3c) are 
shownQ but the subleading contribution of diagram (3e) appears to have been omitted. 
We also disagree with the expression for (3b) in Ref. [^ by a factor of |. However, 
the numerical impact is negligible. 

The corrections involving massless quark loops vanish in the limit mz — ^ and give 
numerically small contributions. For the massless up-type quarks (mi = u and U2 = c) 
we find: 



The leading order term in m|/m? is 



-0.02%|A^,p. 



Sh!g 



N, 



(47r) 



C ,w ,2 



Hjfcl 



rrir. 



rrir 



9m| V K-. 



- h. 



TTlr 



18m? 



(3.10) 



(3.11) 



For niq = 100 GeV, this gives -0.01% lA^J^. 

The massless down-type quark dependent correction is: 



(3.12) 



(3.13) 



For mg = 100 GeV, this gives -0.09% lA^^p. 

Combining everything together, and summing over the generation indices j and k, the 
shift of the i-th generation lepton coupling to the Z due to R-violating interactions is: 



Sh^ = -OM%\K,f 






ig order term in m\/mi is 








' ml ml 
PM8m? ^''^9ml 
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1- 
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-31n^ 



5/^f = E {^^t + ^^ 



ijk 



0.61% E 1^^3.1' - 0.08%E |A:,,,P - 0.08% E I a: 

k k k 

o.6i%Ei^:3.r 



(3.14) 



where we drop the subleading terms. (This is equivalent to keeping only the diagrams 
involving the top and the stop.). 



'c/Eq. 9 of Ref. [10|. 



IV. FITS AND NUMERICAL ANALYSES 

In order to place limits on the A^ -^ couplings from Z decay, we need to know how the 
observables at LEP and SLD will be affected by the shifts bhf in the left-handed coupling 
of Cj to the Z, as well as by other R-conserving vertex and oblique corrections. 

The relevant observables are 

r(Z ^ hadrons) ^ iV^ S,=n,d,s,c.fe(/^L + ^D /, . n 

' r(z -. tn-) {hi + hi) ^ ■ > 



and 



hj -hj 

as well as 

AfbW = \a,A, (4.3) 

where i = e, fi, r. The shift in R^ due to shifts in the coupling constants is 

Re. Thad h^^ + /i£^ 

= A^ + 4.3 5/if (4.4) 

where we have subsumed all the hadronic corrections and the lepton flavor independent 
oblique and vertex corrections into a single parameter Ar, and the coefficient of 5hf is 
calculated for the value sin^ 6^ = 0.2315. 

We note that the A' couplings we are trying to constrain also contribute to A^ through 
Fhad since they modify the Zqq vertices^ as well as the Zii vertices. However, since A^ 
also subsumes highly model dependent corrections from the R-conserving sector, it can be 
considered an independent parameter from 5hf in our fit. 

Similarly, the shift in Ai is given by 

6Ai _ ihi^hjjhi^ -ihlhijhi^ 



A, hi - hi 

- ^A + , 4 ,4 Otlf, 



25 5hf (4.5) 



^The flavor dependence of R~violating corrections to the Zqq vertices can be used to constrain A' 
and A" by looking at purely hadronic observables |11|. 
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where we have subsumed all the lepton flavor independent oblique and vertex corrections 
into a single parameter A^, and the coefficient of 5hf is again calculated for the value 
s\v? 6y[/ = 0.2315. We do not need to introduce another flavor independent parameter for 
Afb(^) since 

^ = ^ + ^ = 2A,-25«,f-25«,f. (46) 

/ipBl^j Ae Ae 

Therefore, we can express all corrections from both R-breaking and R-conserving interac- 
tions in terms of just 5 parameters: Ar, Aa, and 6hf, £ = e, /i,r. A five parameter fit 
cannot be conducted, however, since a change in any one parameter can always be absorbed 
into the other four. We therefore define 



A^^ = A^ + 4.3 ou, , 
A^^ = A^ - 25 6hf, 



Ofie — '^'"^ '^'"e 5 



6re = Shf - 6hf. (4.7) 

and perform a four parameter fit instead. Note that the parameters 

Ik k ) 

{ k k 



,/ |2 



(4i 



are measures of lepton universality violation. The dependence of all the observables we use 
in our fit to the four fit parameters is: 

-^ = Ar.+ 4.3 <5,e 

X ID 

^ = A«^ + 4.3 6re 

^^^ - A.„ - 25 6„ 



A, 



SA 

" A^, - 25 6re 



Ar 

SAFBJe) . 

— j T-- = 2Aa, 

^FB(e) 

-^-^ = 2A..-25 5,. 

^^ = 2A,.-25 5.. (4.9) 
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In table |, we show the most recent data of these observables from Refs. |]T6|, |]T7 



and ||18[. The Standard Model predictions were calculated by ZFITTER v6.21 |T^ with 
standard flag settings for the input values of mj = 174.3 GeV [^, nifj = 300 GeV, and 
dsi'f^z) = 0.120. The limits on lepton universality violation is insensitive to the choice 
of the Higgs mass since Higgs couplings within the Standard Model do not violate lepton 
universality by any appreciable amount. The correlation matrix of the LEP Z lineshape 
data is shown in table ^ 

The result of the four parameter fit to all the data in table | is 

5^e = 0.00038 ± 0.00056 

S^^ = -0.00013 ±0.00061 
Aa, = 0.052 ±0.012 
ArI= 0.0007 ±0.0020 (4.10) 

with the correlation matrix shown in table |T|. The quality of the fit was x^ = 8.3/(12 — 4). 
In figures ^ and ^ we show the la constraints placed on 6^e and 6re in the Aa^ = Ar^ = 
plane by each observable. It is seen that the strongest constraints come from i?^, i?^, and 
Aj. from the r polarization measurement at LEP. In figure |^ we show the 68% and 90% 
confidence contours on the ^^e^f^re plane. 

The limits on (5^e and 6re translate into limits on the R-breaking couplings: 



El^mr-El^ml'= 0.062 ±0.093 

El^m.|'-El^ml' = -0.02 ±0.10 

E l^^r - E l^ml' = -0.083 ± 0.093 (4.11) 

k k 

If we neglect the 13k terms since they are already constrained to be small (recall that 
Eq. p.l3| ) shows the 2a upper bound), we obtain the following la {2a) upper bounds for 



the 23k and 33k terms: 

$:|A'23,P< 0.16 (0.25) 

E|A'33.P< 0.08 (0.18) (4.12) 



or 



IKskl <0.40 (0.50) 

|A:,3,|< 0.28 (0.42) (4.13) 

The limit on A23;, should be interpreted as a limit on A232 since X2^^ and A233 are already 



fairly well constrained by other experiments |T^. If any of the 13/c-terms (in particular, 



A'^32 with a 2a upper bound of 0.28 |jTj]) are non-zero, these limits will be weakened. 

It is interesting to note that since the measured value of Rr is smaller than the measured 
value of Re, this pair prefers a negative value of Sre {cf. Eq. [4.9|). The same can be said of 



the pair Ap^i^) and 74fb(t)- On the other hand, the measured values of A^ (LEP and SLD) 
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and Alr are all larger than the measured values of Ar (LEP and SLD) so these observables 
prefer a positive value of 6re- (This is not apparent in figures | and ^ since they show the 
constraints in the A^^ = plane.) Due to this conflict, the central value of 6re preferred by 
the global fit is virtually zero which satisfies neither the i?'s nor the A's. In fact, A^ from 
LEP, with its smaller fractional error, actually accounts for 2.8 out of 6.8 of the x^ of the 
fit. 

If we perform our fit on the six Z line-shape parameters only, the result is 

5^e = 0.00002 ± 0.00061 

S^^ = -0.00082 ± 0.00070 
Aa, = 0.055 ± 0.033 
Ar^= 0.0022 ±0.0022 (4.14) 

with x^ = 1-9/(6 — 4), and the correlation matrix is shown in table |rV|. The 90% confidence 
contour in the 6^e^5re plane is shown in figures |^ and p. This translates into 





" ) , \^13k — 
k 


0.00 ±0.10 


S^ X' |2 
l^ ^33fcl ' 


k 


-0.13 ±0.11 


ft, 

»r X' |2 

L ^33fcl " 


~z_^W'ik = 


-0.138 ±0.097 



(4.15) 

k k 

Again, neglecting the 13k term, we obtain the following la (2o") limits: 

J2\K3k\'< 0.10(0.20) 

El^ml'< -0.02 (0.09) (4.16) 

k 

or 

l^ml < 0.37 (0.49) 

\K3k\ < (0.30) (4.17) 

The negative central value for 6re leads to a reduced upper bound for | A33^ | . 

If we perform our fit on the two r-polarization observables and the four SLD observables 
only, the result is: 

5^e = 0.0040 ± 0.0046 
5re = 0.0025 ± 0.0013 
Aa, = 0.062 ±0.013 (4.18) 

with x^ = 0.84/(6 — 3), and the correlation matrix is shown in table |V[ The 90% confidence 
contour in the Sfj,e-Sre plane for this case is also shown in figures |^ and |. This translates 
into 

El^ml'-El^ml'= 0.66 ±0.75 



k 
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El^33fer-El^m.|'= 0.41 ±0.22 

E l^mP - E l^mr = -0.25 ± 0.77 (4.19) 



k k 

Neglecting the 13/c term, we obtain the following la (2a) limits: 

ElAkf <l-4 (2.2) 



E|A'33fcr< 0.62 (0.85) (4.20) 

k 



or 



|A'23fc| < 1.2 (1.5) 

|A:,3,|< 0.79 (0.92) (4.21) 

This time, the upper bounds are considerably larger. 

This shows that had we used only the Z line-shape variables, which has been the case in 
previous analyses by other authors [|T0|,12|, or only the leptonic asymmetries, we would have 



reached drastically different conclusions concerning the limits on the R-breaking parameters. 
Only through a global analysis were we able to constrain the parameters in a consistent way. 

V. SUMMARY AND CONCLUSIONS 

We find that flavor- conserwm(7 leptonic Z and W decays can be used to place significant 
constraints on the size of R-parity violating A' couplings. Current bounds from lepton 
universality violation in leptonic Z decays from combined LEP/SLD data are 

|A23fc| <0.40 (0.50) 

|A:,3,|< 0.28 (0.42) (5.1) 

at the la (2a) level, assuming a common squark mass of irtq = 100 GeV and the sup- 
pression of X'l^i. couplings. For larger (common) squark masses the above bounds should 
be interpreted as bounds on |A'| x J F{x)/ F{xq), where F[x) is defined in Eq. ^^ and 

2 

^0 ~ (100GeV)2 • 

Numerically, our numbers are not a significant improvement over those cited in Ref. |]14 
However, the methods used to derive previous limits [050 were intrinsically flawed in that 
(1) R-conserving effects were not properly taken into account and (2) R-breaking effects 
on only the leptonic widths of the Z were considered. Indeed, had we also considered only 
the leptonic widths, our limits would have been those of Eq. ( |4.17| ). The analysis of this 



paper avoids these problems by focussing on lepton universality violation and performing 
a global flt on all LEP/SLD observables. The R-conserving effects are taken into account 
by parametrizing and fltting them to the data also. (Similar methods have been used in 



Ref. [^ to constrain flavor speciflc vertex corrections while taking into account the flavor 
universal oblique corrections.) 

Current bounds on lepton universality in leptonic W decays provides the constraint 

14 



\X',,,\ < 1.7 (2.4) (5.2) 

at the lo" (2cr) level. While not currently competitive with the Z decay bounds, the Fermilab 
results are complementary independent measurements, and they can be expected to improve 
dramatically at Tevatron Run II. 

If the error on the LEP/SLD observables continue to shrink with the current central 
values, then eventually the region allowed by the line-shape variables and the asymmetries 
will fail to overlap in figure ^. In such a situation, not only the SM but the MSSM with 
R-parity violating couplings would be ruled out. In fact, no theory which introduces lepton 
universality violation in only the left-handed couplings would be viable. 

Currently, the LEP and SLD observable provide the best limits on the Agg^i. couplings. 
However, one can potentially place a limit on the A' couplings by looking at invisible decays 
of the T and J/\l/ resonances at the B and r-charm factories |2^. The current bounds on A^33 



imply that the correction to the invisible width of the T resonance can be as large as 30%. 
A rough estimate shows that if the T invisible width is found to agree with the Standard 
Model prediction with 5% accuracy, the A' coupling would be constrained to be |A^33| < 0.16 
at the 2a level. In addition, constraints on IA333I will be available from forthcoming Tevatron 



studies of the decay t — > r6 [23 
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APPENDIX: FEYNMAN INTEGRALS 



Here we make explicit our notation for the scalar and tensor integrals that appear in the 
calculation. The definitions of the integrals are similar to those of Ref. [|^, but there are 
some small alterations made to take advantage of symmetries of the problem. The hat on 
the tensor integrals serves as a reminder of these differences. 



1. Scalar Integrals 



We define the functions Bq and Cq by: 



Co [p^ g^ {p - qf; mi, m2, ms ] = i 



d^k 



{27iY {k^-ml)[{k + py-ml 
d^k 1 



(27r)4 (P _ ^2) [ (^ + p)2 _ ^2 ] [ (^ ^ ^Y _ ^2 ] 



(Al) 
(A2) 



The general form of Bq is given by 

-1 



B,[p';mi,m,\ 



{Any 



m\ ln(m^//i^) - ml ln(m^/^^) , , , ^^ 2 n 

m^ — 1712 



where A^ = a _ j — 1e + ln47r, and |]25 



1 /S 



2 VA 



mt 



F(p2; mi, ms) = 1 + - — - A In ^ - -Vl - 2S + A2 In 



■l-S + V1-2S + A2 



rrin 



with 



2,2 22 

^ ..ii + m^ . m, — m^ 

S = -^-^ — -, A = ^-^5 — -■ 



p 



2 5 



P 



,1-s- V1-2S + A2; 

(A3) 



(A4) 



The function F{p'^; mi, 1712) is well-behaved in the limit p^ -^ 0. For small p^ <^ m\, m^, the 
behavior is: 



F{p^]mi,m2) 



F{p^; m, m) 
F(p2;0,m) 



p^ 



2{m\ — mIY 
^2 



2 2 2m\m\ ( mV 

m, +mo 5 5- In — ^ 

mf — m2 \ mg 



+ 



p 



Qm? 

^2 



P 



2mi^ 



+ 



+ 



(A5) 

(A6) 
(AT) 



Useful special cases of the Bq function are: 
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Bq [p^; m, m] 
Bo[p^;0,m] 



-I 



(47r)2 

-1 

(47r)2 






l_in^^l 



i2' 



Am^ 



-1 



A,-ln— + 2-1- — llnd-Ar 



/x- 



p^ 



m^ 



(A8) 
(A9) 



The general form of the Cq function is fairly complex and we refer the reader to Ref. [^ 
It simplifies considerably for the following cases: 



C'o[0,0,/;m,0,0: 
C'o[0,0,p^;0,mi,m2] 



-^ ^ \\ 


(47r)2p2 _ 

-1 1 


(47r)2 4p2 






Li2 



p^ 



m^ 



\n' 



1 - S - v^- 2S + A2 



In^ 



mj 
mi^ 



(AlO) 
(All) 



where S and A are defined as in Eq. 

The following Cq could be expressed in terms a sum of dilogarithms, but for our pur- 
poses it is simpler to reduce them to a Feynman parameter integral and either perform the 
integration numerically or, if an expansion is needed, to expand the integrand directly and 
then integrate. 



C'o[0,0,/;m,M,0; 



Co[0,0,p^;m,M,M' 



dx 



I 



(47r)2 Jo p'^(^l — x) + rn? 



In 



1 /•! 



dx 



[m? - M'^)x + M^ 
(l-x)(M2-a;p2) 

[rn^ - M'^)x + M'^ 



(47r)2 7o p^{l-x) + {m?-M'^) 



In 



-p'^x{l - x) + M2 



(A12) 
(A13) 



Definition and general form of Bi. 



2. Tensor Integrals 



Ki, 



{2nY (A;2-mf)[(A; + p)2-mi] 



= P^Bi[p'^ -,1711, 1712] (A14) 



Bi [p2; 1711,1712] = --5o [/; mi, ma ] + -—- i ^^ ^ ^ J F(p^; mi, ma) 
^From Eq. (|A5|) , we find that in the limit p'^ -^ 0: 



L^ 



1 11 
5i[0;mi,m2] = --5o [0; mi, ma] + ,. ■.^., ^ 2T 

2 (47r)2 4(mf — ma) 



Other special cases: 



(A15) 



2 2 2rri{ml /m? 
m, + m, ^ ^ In — ^ 

mf — m2 \^2y 



(A16) 



17 



1 1 

(47r)2 2 

p2_»0 1 1 

^ (47r)2 2 



A, 


-In 




+ 2 


9 

nr 


( 

1- 
\ 


p2 


Ae 


-In 




1' 

+ 2 









p^ 



r/i^ 



^i[p';0,m] = ^; A.-ln(^l+2-^-(l-^l In ( 1 - i^ ) (A17) 

(A18) 

(A19) 
= (m^ -ma) -Bo [0; mi, 7712] + {m}^ - mg) _Bo [0;m2,m3] + (mg - m^) i?o [0;m3,mi]. (A20) 

Definition of the C-functions: (Note the difference from the definitions in Ref. ||24|| .) 

(Jb rC 1% ,1 



Useful relations among the 5-functions: 

= 5o[p^;mi,m2] + 5i[/;mi,m2]+Si[/;m2,mi; 



C^[p,g;mi,m2,m3] =i 



(27r)4 (P _ ml) [{k + p)2 - ml] [{k + qf - m^] 



(A21) 



C/,^[]9,g;mi,m2,m3] = i^ 



A-d 



d'^k 



kp^ky 



{2nY (P - mf ) [(A; + pf - ml] [{k + qf - m^] 
= P^iPuC2i + qf,quC22 + (P;,g,. + q^,Pu)C23 + fi'Mi^C'24 

For the purpose of this paper, we will only need to evaluate these functions for p^ = g^ = 
(we neglect final state fermion masses). Q"^ = [p — qY = ~2p ■ q will then be the invariant 
mass squared of the initial vector bosonQ. For this parameter choice, the C-functions can 
be expressed in terms of the 5-functions and C*o as: 



^11 = -7^{^o[0;mi,m2] - BQ[Q'^]m2,m^] - {m\ - m^) Cq} 



Q 

C12 = --^ {^o [ 0; mi, m-i]- Bo[ Q^\ m2, m3 ]-{m{- mj) CqJ 

{d - 2)C'24 = -5i [ Q^; m2, mg ] + {mj - mj) Cn + rn\ Cq 

-g2 (723 + 2(724 = -Bi [Q^;m2,m,] - (m? - m^) (7i2 
= -Bi [ Q"^; m3, m2 ] - (m^ - ml) Cn 

We do not list expressions for C21 nor C22 since we do not use them in this paper. 



(A22) 

(A23) 
(A24) 

(A25) 



We caution the reader that p and q defined here are different from those appearing in the figures. 



TABLES 



Observable 


Measured Value 


ZFITTER Prediction 


Z lineshape variables 


input 
unused 
unused 
20.739 
20.739 
20.786 
0.0152 
0.0152 
0.0152 


raz 

Tz 

<d 

Re 

Rfi 
Rr 

ApBie) 

ApBifJ') 

^FB(r) 


91.1872 ±0.0021 GeV 
2.4944 ± 0.0024 GeV 
41.544 ± 0.037 nb 
20.803 ± 0.049 
20.786 ± 0.033 
20.764 ± 0.045 
0.0145 ±0.0024 
0.0167 ±0.0013 
0.0188 ±0.0017 


T polarization at LEP 


0.1423 

0.1424 


Ae 

Ar 


0.1483 ±0.0051 
0.1424 ± 0.0044 


SLD left-right asymmetries 


0.1423 
0.1423 
0.1423 
0.1424 


Alr 

Ae 
Ar 


0.15108 ±0.00218 

0.1558 ±0.0064 

0.137 ±0.016 

0.142 ±0.016 



TABLE L LEP/SLD observables and their Standard Model predictions. The Z lineshape 
observables are from Ref. ||lH|. The rest of the data is from Ref. |17] and |jl^. The Standard Model 
predictions were calculated using ZFITTER v.6.21 |l|] with mt = 174.3 GeV ^, mn = 300 GeV, 
and as{mz) = 0.120 as input. 





mz Tz 


<d 


Re 


R, 


Rr 


^FB(e) 


^fb(^) 


^FB(r) 


mz 


1.000 -0.008 


-0.050 


0.073 


0.001 


0.002 


-0.015 


0.046 


0.034 


Vz 


1.000 


-0.284 


-0.006 


0.008 


0.000 


-0.002 


0.002 


-0.003 


<.. 




1.000 


0.109 


0.137 


0.100 


0.008 


0.001 


0.007 


Re 






1.000 


0.070 


0.044 


-0.356 


0.023 


0.016 


Rfi 








1.000 


0.072 


0.005 


0.006 


0.004 


Rr 










1.000 


0.003 


-0.003 


0.010 


^FB(e) 












1.000 


-0.026 


-0.020 


^fb(m) 














1.000 


0.045 


Afb{t) 
















1.000 



TABLE II. The correlation of the Z lineshape variables at LEP 
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if^e 


ire 


Aa 


Ar 


^IJ,e 


1.00 


0.53 


0.22 


-0.76 


ire 




1.00 


0.28 


-0.63 


Aa 






1.00 


-0.23 


Ar 








1.00 



TABLE III. The correlation matrix of the fit parameters using all data. 





Spie 


ire 


Aa 


Ar 


ifj.e 


1.00 


0.60 


0.32 


-0.79 


Sre 




1.00 


0.29 


-0.69 


Aa 






1.00 


-0.33 


Ar 








1.00 



TABLE IV. The correlation matrix of the fit parameters using the Z line-shape data only. 





Sfj.e 


ire 


Aa 


^He 


1.00 


0.05 


0.12 


ire 




1.00 


0.41 


Aa 






1.00 



TABLE V. The correlation matrix of the fit parameters using the LEP r-polarization and SLD 
leptonic asymmetries only. 
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FIG. 1. Vertex corrections to W~ -^ Ci^i^i' from R-parity violating interactions. 
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FIG. 2. Wavefunction renormalization corrections to W^ -^ ^Il^i' from R-parity violating 
interactions. 
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FIG. 3. Vertex corrections to Z ^ ^il^il from R-parity violating interactions. 
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FIG. 4. Wavefunction renormalization corrections to Z ^ ^ih^ih from R-parity violating in- 
teractions. 
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FIG. 5. lo" constraints on lepton universality violation. 
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FIG. 6. lo" constraints on lepton universality violation. Note the larger scale with respect to 
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FIG. 7. Confidence contours from all data (68% and 90%, in gray), the Z line-sliape data only 
(90%, dashed line), and the asymmetry data only (90%, dot-dashed line). 
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